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REPORT No. 21. 

THEORY OP AN AIRPLANE ENCOUNTERING GUSTS, II. 1 



By Edwin BmwBiL Wilson. 
INTRODUCTION. 

This discussion is an immediate continuation of my previous 
treatment of the subject published in the First Annual Report of 
the National Advisory Committee for Aeronautics, Washington, 
1915, pages 52-75 (S. Doc. 268, 64th Cong., 1st sess., reference to 
which will be by pages). The notations of that work will be con- 
tinued without change except as hereafter noted. 

PERIODIC LONGITUDINAL GUSTS. 

That there is a certain degree of periodicity in gusts is obvious 
from casual observation, from the records of scientific observatories 
like Blue Hill, and from the familiar fact that- all such phenomena 
in nature reveal a general tendency toward periodicity. Needless to 
say the periodicity is not mathematically exact in its regularity nor 
indefinite in con tinuan ce. 

The object, however, of an investigation of the effect of periodic 
gusts on an airplane can for practical purposes be no other than 
to reveal any exceptional effects that periodic, as compared with 
single, gusts may have upon the flight of the machine; and these 
exceptional effects will probably be indicated with sufficient practical 
completeness by an analysis built on the assumption of strict peri- 
odicity, long continued in operation — the phenomenon most to be 
feared being resonance. 

The longitudinal gusts are in, 1°, head-on velocity t^; 2°, vertical 
velocity w t ; 3°, rotary velocity q t . Very little is known as to the 
nature of rotary gusts (p. 65) and hence 3° may be left aside. It is 
not easy to see how vertical gusts can have any pronounced peri- 
odicity; the disturbance of the airplane's motion by vertical gusts 
is (p. 64), except for very sharp gusts, essentially a convection of 
the machine with and by the gust; for both these reasons 2° may be 
discarded. This leaves only 1° — periodicity in the head-on gusti- 
ness — as likely to be of interest. 



i Reprinted from Proceedings American Philosophical Society, voL M, 1817. 
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The gust may be assumed in the form 

J sin pt or u^JeW. (1) 

The differential equations are (p. 59) 

f(D)u = -(0.128Z?»+1.160Z) 1 -{-3.385Z>+0.917)u 1 , 
f(D)w*=-D 3 (0.557D+2A58)u u (2) 
f(D)d = -0.02851Z?m m 
with/(Z>) =Z>* +8.49Z?" +24.5Z? 1 +3.3S5Z? +0.917 

= (Z>*+8.359Z> +23.37) (7J , +0.1308Z> + 0.03924). 

In the previous investigation it was found that the short-period 
heavily damped oscillation was not of much significance except in 
the case of a sharp up-gust (pp. 62-69), and that its significance in 
that case was not revealed in the major motion of the machine but 
in the initial acceleration (or stress) upon it. It may therefore bo 
expected that for periodic head-on gusta the short-period motion 
will be negligible in its effects. It is consequently desirable to carry 
out the numerical analysis in such a way as to separate, so far as 
may be, the short and long natural periods of the machine. 

Let us separate into partial fractions the operator 

I i 
/(2J) = (D 2 +8.359 D +23.37) (D 1 +0.1308D +0.03924> " 



or 



I 0.016Z? +0.089 _- 03JI601Z? +0.24263 
j(D) "* IP +8.359Z? + 23.37 + LP +0.1308U +0.03924" 



(3) 



The first fraction has to do with the short, the second with the long 
oscillation. The two operators are to be applied to certain expres- 
sions derived from (i) by substitution in (2). 

If D=ip, the numerators of (3) have the respective magnitudes 

(0.089* +0.016 , p , ) I /» and (0.0426» +0.016 , 2>')V. 

For p=0, the second is about half the first; for p = °°> the two are 
equal; the numerators therefore do not differ greatly in magnitude 
for any value of p. 

The ratio of the denominators is 



„ T (.03924-g*)' +.1308"?' I 1 /* 
• K = L (23.37 -p»)» +8.329^ J ' 



and is very small when p is less than 1. For larger values of p, we 
have approximately 

l/iP = 1 +23.2/P* +545/p 4 . 

Hence the short oscillations may be neglected when p < 1 without 
introducing much error; but as p increases beyond the value 1, the 
importance of the short oscillation grows rapidly. 
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Consider first the case p<l, neglecting the short oscillation. 
The particular solutions for u, w, and fl, that is, I u , I v , I a , are ob- 
tained from the imaginary parts of 

0 -.01601pi +.04263 , noo „ 
/ = + .1308gi+.03924 ( ~- 0285 

To estimate the value of v corresponding to the maximum dis- 
turbance we may examine the amplitude oid/J, which is 

e nOQKtm T (.04263)' +(.016p)» 7/» 
amp. 7 - .02851J) [_ CflgflaW+OM o8p). J ■ (5) 

The calculation gives 2>»=0.0394 or £=0.1985. The value of the 
amplitude is then about 0.0095/ radians or 0.54/ degrees. If J 
should be 20 ft./sec, the forced oscillation would have an amplitude 
of about 10°. 

As the use of p =0.1985 in calculating is somewhat more compli- 
cated than the use of p=0.2, and as the change from 0.1985 to 0.2 
does not materially alter the amplitude of the forced oscillation 
(and probably does not exceed the error of observations), we may 
use p = 0.2 in calculating the effect of a periodic gust of maximu m 
resonance on the airplane. We shall first note that for y=0.2 
the ratio of the amplitudes of the two fractions in (3) is of the order 
400 to 1 ; and the first fraction is therefore entirety negligible in 
determining the particular integrals. 

For the seoond fraction we have" the complex value 

4.263 - .32* (4.275,- 4.3°) Z4.275, », n o\ 
2.616i - .076 = (2.617,91.6°) "A2.617' y5,y /' 

where the parentheses contain the polar coordinates of the complex 
numbers. The expressions into which this is multiplied to determine 
the coefficients of e { ** are for u/J, w/J, 6/ J, respectively, 

-0.922-0.676i= (1.144, 216.24°), 

0.0983 +0.00456i= (.0984, 2.67°), 
- 0.0057* = (.0057, - 90°) . 

Hence the values of u/J, w/J, 6/ J are 

«//=(- .965 + 1.65i) (cos M +i sin J&), 
w/J = ( - .00918 - .164i) (cos j8«+isin^), 
fl/J = ( — .00948 + .00098*) (cos £t+ isrin 
and I U =J (1.65 cos J&-.965 sin M), 
I V =J (-.164 cos Jgf-.0092 sin JBt). 
I, =>J (.00098 cos .#-.00948 sin M), 
I f '~J (-.0019 cos J#-.0002sin JBi), 
7 W =1.65.7, J„=-.164/, = . 00098/, Z'i,= -.0019/. 
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On substituting these values to find the constants of integra- 
tion (p. 61), it is found that A and 0, corresponding to the short oscil- 
lation in u, are negligible. Also S= — 1.65J, D = .72&J. Hence 

u = Je-«»«(- 1.65 oos. 187i+.726 sin .1870 

+.7(1.65 cos £t-M5 sin ,2t). 

In like manner (p. 62), A' and C are small and B r -.176/, D' 
= -.051./'. 

w=J(r-°^*(.l76 cos. .187<-.051 sin .1870 

-J(.164 cos .£f+.009 sin M) - .0l2j<r*-™ cos 2.43«. 

(The last term is added as a oheck on the initial condition w = 0.) 
Finally (p. 62), ^" = .00007/, £" = .00104/, i7" = .0109J, and 

e=Je--<*s*t (-.00104 cos ,187i+.0109 sin .1870 

+J (.00098 oos 00948 sin M) + .00007 Je- ,JW cos 2.43f. 

Now to find the rise of the machine when the gust strikes it (p. 64). 

w + 1 15.59 = Je- 06 "* (.056 cos. 187* + 1.208 sin .1870 

-J (.051 cos .£« + 1.064 sin .St). 

The cosine terms may be omitted. The integration then gives 
2 = 5.32J cos M+0.44J-J<r M4t (2 sin .1872+5.76 cos .1870. 
A table of values of z may be computed as: 

* = 0, 2, 4, 6, 8, 10, 12, 14, 

2/J=0, 0, -.15, -.54, -1.16, -1.90, -2.60, -2.97. 

This shows the rise or drop, according as J is negative or posi- 
tive, during the first quarter minute. The values of z now fall off, 
pass through 0, and only become large as f nears 35. The natural 
osoiUation is then becoming less effective relative to the forced 
oscillation which has a double amplitude of about 10.6J, or 202 ft. 
if J =20 ft./sec. 

As the existence of a regular periodic gust for any long time is 
almost unbelievable, the only real interest in the calculation is in 
showing that- during the first 15 seconds the effect of resonance fails 
to become so far established that the motion differs appreciably from 
that due to the simple head-on gust previously studied (p. 74). 

In the case of the machine constrained to remain horizontal during 
flight (by some automatic steering device), the corresponding equa- 
tions (p. 69) are for u 1 =Je i ^ t 

u_ .128j» + .598 v 
J~~* .598 -p*+4.078pi 3 ' 

% ^Z£» Le<pt 

7 .598-p* + 4.078pi ' 

As the natural motion is no longer periodic, there can hardly be 
any such thing as resonance, in the usual acceptation of that term. 
We can, however, ask what value of p will make wjj a maximum 
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and. hence induce the maximum oscillation in the vertical motion. 
To maximize 

(MS-p^f+Am&p* ° r 4.078* +(p-.598/3>) 1 
take p* = 0.598 or p = 0.774. The value of w/J is then 

w/J= -0.136e*J* 

and the amplitade of w is 0.136J. The amplitude of the oscillation 
corresponding to the particular solution is 0.1 75J. 

Thus again it is seen that the steering device makes the motion 
far easier than when the machine is free (p. 70). There seems to 
be no need of carrying out the details of the integration. 

NOTE ON RESONANCE. 

_ In defining, by implication, a state of resonance in the calcula- 
tions above, I have assumed that it was the angle 6 which was to be 
maximized by the proper choice of the frequency p of the applied 
periodic force. It may be well to take up the theory of resonance in 
a little greater detail, for there are complications in the kind of 
system we have here to consider. 

A. G. Webster, in his "Dynamics of Particles, etc.," Teubner, 
1904, page 175, gives general formulas for resonance and shows that 
if the damping coefficients are small and if the frequency of the 
impressed force nearly coincides with that of the natural oscillation, 
the amplitude of the forced vibration, will be relatively large. 

_ This is not enough. For in the first place the damping coeffi- 
cients in the case of the aeroplane can. hardly be regarded as small 
(they sometimes exceed the frequencies); in the second place, we are 
not even certain that the motion of the system is wholly oscillatory 
(some of the roots may be real, and even positive if the machine has 
a certain amount of dynamical instability) ; and in the third place, 
under such conditions, the amplitude of the forced oscillation may 
be considerably greater when the frequency of the applied force is 
materially different from that of the system (supposed oscillatory) 
than when the system and the force are nearly synchronous. 
The ordinary theory of simple resonance depends on the equation 

(,D t +'kD+n)x=Jsmpt. 

The particular solution 

is the imaginary part of the expression 

n—p'+ipi' 

The amplitude of I f is the same as the modulus of the complex 
value x. The modulus of is 1; that of z is 

amp - ^[(n-tf'+fcpf"' 
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To make the denominator a minimum we have merely to minimize 

(n-j^+Pjj, 2=p»>0. 

We find. (7=7i— JP, necessitating n>ik t . If, then, the maxi- 

mum amplitude of 7, is 

max. amp. 

where the positive or negative sign must be taken according as Tc is 
positive or negative. It the maximum amplitude for I x 

occurs when p = 0 and is J/n. 

The amplitude is large when Tc or (n— Ji 1 ) is small; it is very 
large when both conditions are satisfied. The largest possible value 
occurs when n=\Tc* and is -fij^c*. In this case the applied force 
has an indefinitely small frequency where the natural oscillation has 
the frequency £/V2. The theory of the system here considered is 
given by Webster (op. cit., p. 155). 

The case which corresponds to that in which we are interested is 
where the system starts from rest at the position of equilibrium. 

The motion is then denned by the equation 



Under normal conditions this quantity remains tolerably small 
until the natural motion is nearly damped out or until that motion 
has time to increase greatly 0c > o). Even if n= iJP+flc*, the 
equation becomes 

a^^rCe-*" cos ifa-cos €&f)+|r (*-** sin \U - sin Jet), 

and the conclusion still holds. 

For the simple system ordinarily treated for resonance the state- 
ment that the motion must be only slightly damped and the frequen- 
cies of the natural and forced vibrations must be reasonably near 
together, is therefore amply justified. The result holds even when 
n < \Tc*, in which case the maximum amplitude for I x (resonance) 
occurs when p=o and is J (n. 

The next simplest case is like that which arises in treating the 
constrained longitudinal motion (0=o) of the aerpolane (p. 69): 

(D+a) u+lw~ — a«j — tw lt o=.128, I— — .162, 
cu+(D + d)w = — cttj— dw t , c=.557, <2=3.95. 

The natural motion is given by 

A'=D>+ (a+d)D+ (ad-bc) =o, 
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and in this case by IP +4.078.Z? + .598 =o. Here the roots are 
both real, viz.,— 3.93 and —0.15. So far as the equation in D is 
concerned we have the case where h is large and n is small. Hie 
equations for the forced motion are 

A'u= — (aD+n) Uj^—lDw^ 
A'w = — (dD +n)w l — cm v 

The question now arises: What is it that is to be a maximum? 
For some purposes it might be the variables u or w — for example, 
the whole theory of gusts here given depends on the gust being 
small and producing small effects, and if by an applied force, the 
values of u or id should become too large, the theory would become 
worthless. Again, if the question had to do with the strain on the 
machine, the derivatives duldt and dw/dt would be the essential 
objects of interest, and should be maximized. Finally it might be 
the values x=fudt and z=fwdt — the actual displacements of the 
machine — which we desired to examine. Let us therefore consider 
several problems seriatim. 

Case 1. — To maximize u with a head-on gust u^e^. 

u ofr+n y .128*? + .598 M 

U A' € *~ .598 -t +4.098*/ * 

The Trm-rimum value of 

.128 a p»+.598 8 .128 a (p* +21.83) 
^9lf = p*F+^098 1 p i 2>«+15.59p J + .3576 

occurs when p* is o, that is, "resonance" occurs for p=o, the ampli- 
tude of the force and the oscillation being the same. 
Case S. — To niaximize w with a headgust. 

This was treated above (sec. 10). The ratio .136 was found; the 
required value of p was .776. 

Vase 8. — To maximize u with an up-gust w v 

.162pt 

u " .598 -p* +4.098?** ' 

The condition is p=.77Q as is case 2; the ratio is .04. 
Case 4. — To maximize w with an up-gust. 

3.95pt + .598 „ fBt 
.598-p*+4.098pi " 

The ma-n'TmiTTi value of 

3.95 g p* + .598 a 3.95»(p*+.Q228) 
(.598 -p*) 2 + 4.098 J p* "p* + 15.59J? 2 + .3576 

occurs when = .022 and p = .15, and the amplitude ratio is about 1. 

We note the very different values of p thus found, namely, 0, 0.15, 
0.776, according to the choice of case. If in case 1, we had taken 
p=.15, the amplitude ratio would have been about .7 instead of 1* 
if p=.776 had been assumed, the ratio would have been .37 instead 
of 1. 
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Case 6. — If we desired to maximize z we should have had to treat 

1 3.95i+.598/j? w 
? .598 -p* +4.098^ ' 

which would have given an infinite amplitude ratio for p = 0. 

Now if we turn to the free machine and try to maximize fddt 
instead of 6, we have to maximize 

.04263' + .016»j>* 
(.03924 1 -p? + .1308 , 2)* 

instead of (5, sec. 6). The value of p* is about .0307 and of p about 
.175 instead of .2 as before. The amplitude ratio is then only 
slightly in excess (about 4- per cent) of that previously found — in 
other words the numerical values are such that resonance for 0 and 
for fddt, which is the preponderating term in the expression for z, 
occurs for considerably different values of p, but the effect is about 
the same. This may be regarded as validating our procedure 
(sec. 6) in maximizing 6 instead of fddt 

To sum up this discussion of resonance as applied to the aeroplane 
we may say that the frequencies which produce "resonance" depend 
largely upon the quantity in which the effect of resonance is to be 
sought and that the frequency which makes for a strong resonant 
effect in one quantity may make on another an effect much weaker 
than the maximum — or it may not. 

There remains to discuss the question whether the effect of reso- 
nance is practically serious, i. e., whether as in the case of the motion 
of the machine, above treated, the effect fails to make itself felt 
until after so long a time that the pilot would be entirely able to 
deal with it or the wind would really have in all probability ceased 
to be periodic with the period required. 

Now in order to insure that resonance is effective both of itself 
and as against the natural motion, we should reasonably oxpcct to 
require, 1°, that the resonant frequency p be large (for if it be small 
the pilot will have ample time to take care of it), and that, 2°, itrbe 
reasonably different from any natural frequency which is only 
slightly damped (for in the latter case the initial conditions will 
probably be such as to cause the natural and forced effects to clash 
for a considerable interval of time). 

This problem in its generality is so complicated that I have as yet 
been unable to determine whether there may be practically serious 
effects due to resonance, but from the cases I nave here treated, from 
the general considerations which I have advanced, with duo regard 
to the restrictions on p which appear to be reasonable, and from 
cases which I have examined without mentioning them here, I 
should judge that resonance is not a practically serious matter in 
longitudinal motion, and that we may safely confine our attention 
to gusts of the form J" (1 — er*). 

One type of resonance which deserves consideration is that of the 
damped harmonic gust JtT* sin pi. It would be conjectured that 
if — « ±pi were nearly equal to a pair of roots of A = o, there might 
arise a considerable disturbance. It is not likely that a gust of this 
type would exist in reality, but the commencement of any gust 
might resemble very closely the commencement of such a gust and 
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if tiie effect of this type -were very marked as compared to that of 
the types already considered, it -would be necessary, for the sake of 
foreseeing the worst that could happen, to discuss this type. _ 

I have not time to take the matter up here. Moreover, I imagine 
that it would be found that the constants of integration turned out 
to have such values that the gust, though tuned in damping and in 
frequency to the natural motion of the machine, did not have very 
large effects except in cases where n and p were small enough to 
allow the pilot easily to correct for the disturbance. 

The damped periodic gust has been treated by Brodetsky, 1 who 
finds the amplitude of the particular solution is a maximum (for the 
machine I am dealing ■with) when t= 16 sec. and is thena tolerably 
large quantity, — but the pilot has a quarter of a minute in which to 
react to his environment. It is, however, by no means certain that 
the pilot would have to react so quickly — the constants of integra- 
tion might turn out, as I have fust suggested, such that the motion 
during the first quarter minute was not far different from that in 
the case of the simple gust. This was what was found to happen in 
the case of the periodic gust above treated (sec. 9). The amplitude 
of the vertical motion so far as the particular solution was concerned 
turned out to be about 5.3 J, but the constants of integration were 
such as to postpone the major effect of the particular solution until 
30 or 40 seconds had elapsed. If we have a damped harmonic gust 
and such a postponement were operative, the damping would become 
effective and the gust might turn out to have at no time an effect 
much in excess of the maximum effect of a single gust of the form 

/(l— «-**). 

INFINITELY SHARP GUSTS. 

In my previous paper I discussed gusts «T(i — e - ") rising from zero 
to J with various degrees of sharpness depending on the value of 
r — the larger r, the sharper the gust. An mfinitely sharp gust would 
be one for which r was indefinitely large. Such a gust would repre- 
sent an absolute discontinuity in the velocity of the wind. This is 
impossible, though it represents a state of aerial motion which may be 
nearly approached. Moreover, the ifinitely sharp gust could not 
strike the machine all over at once, and hence the theoretical effect 
of such a gust on the assumption that the machine is instantaneously 
immersed must differ from the actual effect upon a machine running 
into a discontinuity in the wind velocity. 

For this reason one may well limit his considerations to finite gusts 
with a value of r not greater than 5, say, as I did. Nevertheless 
if the calculation of the effect of an infinitely sharp gust is simpler 
than for a finite gust and if the limiting motion derived for such a 
gust is not appreciably different from that for a sharp gust of rea- 
sonable sharpness, the discussion of the limiting case will be justified. 

Consider first the longitudinal motion and a head-on gust ii* 1 =«7'- 
(i— e"*'), r enormously large. According to the symbolic method 
D= —r must be substituted to find the particular solution for « -rt . 
As, however, A is of the fourth degree in I? and all the polynomials 
upon the right hand are of degree 3 or less, the result of the substitu- 
tion is easy to find. 



t ASrcnautfcal Journal, London, 20, 1915, p. 1H. 
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Eor example, when u±=*J(i— e - "), 

IJJ--I- <r^.(.128/r), 1^-= -J, 

/„//= - <r^(.557/r), /„ =o, 

Id/ J = - e^(.02851/r») -0, 70 o = 0, 

7' fl /J"=ff- rt (.02861/r I )=0, 7'<? c =0. 

The equations of motion are 

u/J=, g-*- ts %0009 coa 2.43* +.0032 sin 2.43*) 

+ «-■<>«« (.9991 ^ .ig7* + .3577 sin .187*) -i -«-*(. 128/r) 

u,/jr =< r<.iM(,1066 cos 2.43* -.0435 sin 2.43*) 

+<r-<««(-.1066 cog .i87i + .0352 sin .187*) -dr"(.557/r), 

1000/J (-.0402 cos 2.43*-.0278 sin 2.43*) 

+«r-«««(.04O2 cos .187*-.6683 sin .187*). 

The calculation of the constants of integration is much simplified. 
The terms er^lr are retained because the stresses (forces) due to the 
gust are calculated from du/dt and dwj&t to which these terms make 
an initial contribution — there is an instantaneous initial stress. 
When* = 0, 

du/dt - (. 128 - .004 - .008 - .085 - .067) J = - .016/, 

dw/dt = (.557 - .446 - . 106 + .007 + .006) J = .01&7. 

These are the initial accelerations and should vanish because the 
gust though infinitely sharp begins at zero. That they do not 
vanish is due to an accumulation of errors. 

Immediately after the initial instant, however, the first terms, 
viz., .128 and .557, being multiplied by e^ 1 vanish. The other terms, 
however, being multiplied by comparatively slow changing functions 
are not altered. Hence immediately after the first instant there are 
accelerations — .128J and —.557 J along the x and z axes, respectively. 

To put it another way, there is a discontinuity in the stress at the 
initial instant — as might be expected. The amounts of the discon- 
tinuities are also just what might be expected, viz ? X U J and Z V J. 
In like manner for an up-gust the initial discontinuities in accelera- 
tion are X W J and Z U J along the % and z axes. These results could 
have been foreseen from the differential equations themselves as 
well as from "common sense." The path in space is not materially 
different for an infinitely sharp gust from what it is for a reason- 
ably sharp gust. ' 

It-may therefore be said that a tolerably good idea of what hap- 
pens for sharp gusts may be had from the consideration of infinitely 
sharp gusts. 

It has just been stated that the conclusions concerning the initial 
accelerations may be foreseen from the differential equations. This 
may be proved as follows: We have 

(7? - X u )u - X a w - Xm -g3 = XJ& + X^ + X$ lt 

- Z u u + (D-Z a )w-(Z q + V)q = Zji x + Z v w y + Za u (6) 

- M u u - M u w + Qc B 2 D - M q )q = + + JLg,, 

770-2=0, 
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where the equations have been reduced to four involving only the 
first derivatives of the four variables tt, w, q, 6, with the initial condi- 
tions tt=u?=<7=0=O, by the device of choosing <[=DQ as an inde- 
pendent variable so as to eliminate the second derivatives. 

These equations determine the first derivatives at the initial 
instant or at any instant in terms of the values of the variables at 
that instant, namely, 



(7) 



At the initial instant u, w, g, 6 vanish. 

"With an infinitely sharp gust u u w u q t may be considered as not 
vanishing, but as starting at finite values, J„ «/„, J q . The derivatives 
are then at the initial instant 

Dw=Z u J a +Z u J a +Zjg, (8) 
TcjDq- jq,J u + M a J a + MfJa, 
2)0=0. 

The first two equations give the X and Z accelerations of the 
machine which determine the stresses as the accelerations times the 
mass. 

We have, for numerical values, 

Du= - .128J a +.162J„+0j;, ifX ? = 0, 
Dw= - .557J v -3.Q5j a +QJ a , if Z„=0, 
342? 2 = 0 J M + 1.747,, - 150 J q , if 2£ = 0. 

The last equation determines the couple. tending to break the ma- 
chine, by bending in the »-2-plane, on multiplication by the mass m. 

That which I hare called an infinitely sharp gust is not an impulsive 
gust. The impulsive gust is both infinitely sharp and infinitely 
intense, but endures for only an infinitesimal time. The effect of 
an impulsive gust is to produce instantaneous changes in u, w, 
Such an impulse, like the impulses of ordinary mechanics, puts an 
infinite strain on the machine for an infinitesimal time, and the 
only way to tell whether the machine will stand the strain is to take 
the yielding of the framework into account — it is a problem in elas- 
ticity. For the purpose of calculating the stresses produced by gusts 
on the machine l therefore prefer the sharp gust to the impulsive 
gust. 

For the purpose of treating the motion of the machine after the 
gust strikes it — the gust being now a sudden fierce squall in other- 
wise still ah* — we have merely to determine the constants of integra- 
tion from the initial condition %» a? 0 , g 0 , and 0=0, where u 0 , w 0 , q 0 
are the impulsively generated velocities. These equations are (p. 61) : 

u 0 =A+B, 

w a = -4.04A+34.5C-. 10582? + .0025872?, rcn 
0 = - .132A- .0946(7+ .0024782?+ .0057992?, w 
2o = .7034 + .205<7- .0012462? + .0000842?. 
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Analytically the effect of the impulsive gust upon the equations 
for determining the constants of integration is merely to replace the 
initial values of the particular solutions I m I m , 1'^, obtained 
on the hypothesis of finite gust9, by the respective values — u^, 
— w 0 , 0, —2,. The effect of the disturbance may therefore be calcu- 
lated at once from my equations (23), (24), (25), (26), as soon as the 
values u 0J w oi 2q have been determined. 

In the calculation of u„, w 0f q t , the same doubt arises as in the 
theory of any very sharp gust, namely, the effect of the partial 
immersion of the machine, is the effect of a blow traveling along a 
mechanism the same as that of the blow applied instantaneously at 
all points of the mechanism? The possibility of a difference be- 
tween the instantaneous immersion and the immersion distributed in 
time would arise only if, 1°, the machine had time enough to change 
its orientation appreciably or, 2°, the acquired velocities were suffi- 
cient to change the relative wind and thus affect considerably the 
impulsive pressure. 

Even if we assume that no material difference in effect is to be 
expected, it-is difficult, to make the proper assumptions to lead to 
reasonably satisfactory values for «,>, w 0 , ?» for any actual machine 
whose characteristics are expressed in terms of tlie mechanical 
coefficients m, fa*, Z7, and the aerodynamical coefficients X u , X w , X v 
Z„, Z a , Z q , M u , M w , M q . It is by nojmeans certain that for a "con- 
siderable aerial disturbance the finite instantaneous changes in u, w, 
q can be calculated from the equations (8) by replacing D by the 
sign A for the increment and taking J U} J„, Jg as the intensities of 
the impulsive gusts; for the nine coefficients X M , etc., vary with the 
intensity of the relative wind. 

It is for this reason that I have used finite gusts of various degrees 
of sharpness instead of impulsive gusts. Moreover, it is not certain 
but the finite gust represents more nearly actual conditions in the 
air when flying is at all possible. 

An article by Brodefcsky, with an introduction by Bryan, has re- 
cently reached this country, 1 in which impulsive gusts are considered, 
relative to Bryan's skeleton aeroplane consisting of a forward main 
plane and rear tail plane. The discussion is both interesting and 
important as is everything to which Bryan, the great pioneer in this 
subject, sets his name, but it does not seem to help me, so far as I 
have yet been able to examine it, in regard to the effect of an im- 
pulsive gust upon a machine whose properties are actually de- 
termined in the wind tunnel. I have therefore decided to let stand 
the brief general considerations above. 

THE ACTION OP THE AIR SCREW. 

In the work to this point, I have made for the discussion of gusts 
the same assumption concerning the action of the propeller that 
Hunsaker, Bairstow, and others have made for discussions of stability, 
namely, that under varying conditions the motor speeds up or slows 
down so as to deliver a constant thrust along the af-axis. 

It would be equally reasonable, from some points of view more 
reasonable, to assume that under changing conditions of relative air 
velocity a motor speeds up or slows down so as to deliver the same 



i Aeronautical Journal, London, SO, ISIS, 139-15*. 
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effective horsepower. We should then have the power P equal to 
the thrust 22* (taken positive) multiplied by the velocity — U: 

P= -2277= - (H+dH)(U+w), 
UdH+uH=Q, 

dH—H^—P-^. (10) 

This is an additional force which is directed along the X-axis if 
the propeller shaft is horizontal for the velocity of flight — U. If 
in the standard condition the shaft is not horizontal there would be 
components 

-P pjcos a, +P-p£in a 

along the x and z axes, a being the angle from the horizontal up to 
the direction of the shaft. Purthermore if the shaft did not pass 
through the center of gravity there would be a pitching moment 
—Phu/TP if h is the distance of the line of the shaft above the 
center of gravity. 
The equations for the natural longitudinal motion would then be 

(i>-x.+^5)tf-Zww-cav?+rfff-o, (id 

the other two equations remaining unchanged, if we assume for 
simplicitv that a=7&=o. The effect of the varying thrust is to 
ohange X« to X^—Pg/mTP. We have the value .128 for 

this machine. If the effective propeller horsepower were 87 for 
Z7= - 115.5, the value Pgfm TP is 

Pg 87X550X32 
mt> _ 1800 X 13350 = ' 

The modification of the equations of motion on replacing X u 
= —.128 by X a = — .191 would make an appreciable, though not 
very serious ohange. 1 
The determinant A would become 

342?<+290.82? s +850.92? 1 + 165.l2?+31.18 

=34(2?* + 8.5532?'' +25.032? 2 +4.8562? + .917) 

as compared with 

34(1?* + 8.490.0 s +24.502?* +3. 3852? +.917). 
The rapidly damped oscillation would, as a first approximation, be 

-4.276 ±2.596i instead of-4.245 ±2.545*. 

The first approximation for the small root would be 

-.097±.177i instead of -.069±.18l£. 

The damping would be more pronounced and the oscillation a trifle 
faster. 

It may be concluded that whether the screw is supposed to 
deliver a constant thrust or a constant power is not very important 
to the theory either of stability or of gusts. It is not unlikely that 

29165°— S. Doc 128, 65-2 27 
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the actual behavior of the screw lies within the limits set by these 
two assumptions or sufficiently near to one of the limits to validate 
the use of either hypothesis. 

The Aeronautical Journal; London, 20, 1916, p. 142, quotes 
Bairstow and Fage as giving the formula 

dH= — .OllHdV, ' V in miles per hour, 
which is dE= — .0073fli£F, V in feet per second. 

With 27= 115.5 numerically we would have for constant power ■ 

dH=- .0086653 V, V in feet per second ; 

and, if I understand correctly the use of the signs + and — in the 
quotation, the results are in as good agreement as could be expected 
in view of the fact that I have no Knowledge of the value of JJ 
for which the data quoted are given. (If the motor and screw were 
exactly designed to give a maximum efficiency at a standard speed 
U, we couldnot expect the efficiency to be the same at relative air 
speeds either higher or lower, and this would slightly influence the 
result.) 

EQUATIONS FOB LATERAL MOTION. 

The differential equations for the lateral motion of a machine in a 
gust may be written as (p. 54): 

dv/dt+g<p+ Ur-Y&- Y p p- Y>= Y,v t + T^+Y^ u 

A/m. dp/dt—L„v— Lpp — = L,v t + L p p l + Lf^ , (12) 
C/m. dr/dt — N,v — N p p - Njr = N v v t + N p p x + 

where the terms mvolving the small unknown product of inertia E 
have been neglected and gusts of the type v lt p 1} r x have been 
allowed. 

The gust v K corresponds to a side wind. A change in the direc- 
tion of the wind by a small angle would produce such a gust even 
in absence of any change in the wind velocity. The gust p, is a 
rotary gust tending to produce a bank: as a disturbance in the air 
itrwould correspond to a horizontal roller run into end-on (axially). 
The gust r t corresponds to a column of air rotating about a vertical 
line. 

This last is a common type of aerial disturbance, easily observed 
on a warm day, often of very small diameter compared with the 
spread of the wings of an aeroplane, and accompanied by a strong 
rising current of air. Such a vertical vortex, if small, might strike 
one wing of the machine alone, and, due to the rising current, heel 
it over suddenly. It is, however, not this small local disturbance 
which we can consider by our methods here, but the larger and more 
gentle rotation in the an* which might immerse the whole machine 
many times over and which produces a yawing motion in the 
machine rather than (primarily) a roll or bank. 

31. Place D==d/dt. Then the equations are 

(Z>- ?>+ to- YpD)<j>+(U— F r )r= Y&+ Ya+ 3> u 

-X,v+ (kJD-L p )l4-Lr-L&+L^+£r u (13) 
- N,v- NpD*+ (Jc/D-N r )r= N 9 v t + N^+N^, 
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where Jcs*-=A/m and Te<?=0!m. The determinant whose vanish- 
ing determines the natural motion is 

D-Y, ff—YpD 
-L 9 W-LpD 



A = 



17- Fr 

-A- 

7c<?D-N T 



Let the cof actors of A he 
-NJ) Tc<?D-N r 

-Lr 

1cSD-N r -«„ 

-N, -NpU 
-NJ) lc<?D — N r 



5„ = 



8„ = 



* 259223* + 2314023 3 + 847823. 



-£« | = -59.5523^26.55, 



g-Y p D U-Y r 



D-Y, U-Y r 
-N, Te<?D-N r 
. I -tf, -^23 

^-FJ? U-T, 



= -32.842? 3 - 280.723, 

227023-868.8, 

| =70.6Z3 S +44.5D + 109.9, 
=28.76, 



5» = 



. 4243Z3* + 36270Z3 - 1776, 



.36.723 s +323.lZ3 1 + 77.88Z3 + 27.15, 



IP—LpD -1^' 
Y.j^ T ^ZxJ 1 \ =55.223+ 111.2, 

where the numerical values are those arising from the data de- 
termined for the Curtiss Tractor (which is tine maffhnm under in- 
vestigation) by Dr. J. 0. Hunsaker as given on page 78 of his paper, 
"Dynamical Stability of Aeroplanes, Smithsonian Misc. Collect., 
Washington, Vol 62, No. 5, pp. 1-78, 1916, namely, 

Y p =0, Y T =0, 
L p = -314, Iy= +55.2, 

4=0, tf r =-27.0, 

fc<,'=70.6-, Z7= -116.5, y=32.17. 

The value of A is then (23- Y,)d 11 +gd i2 + US U or 

A = 259223*+ 2378023 8 + 1800023 s +3461023 - 854. 

This result checks with Hunsaker's doc. cit., p. 78) as well as prob- 
able. The equation A = 0 may be written as 

23*+ 9.17223 s + 6.94323* + 13.3523 - 0.3295 = 0. 



F,= -0.248, 
0.844, 
N,= -0.894, 

^=36.7+ . 



32. From the last two terms, one root is indicated 
=0.02468; and the correction can readily be found, giving 

23=0.02436. 

There is another root near 23= —8.5, the exact value being 

23= -8.542. 



23 
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The other factor of the biquadratic equation is 
D l + 0.65370 + 1.583=0, 
of which the roots are 

0 = - 0.326.8 ±1.215i. 

The complementary functions for v, 4>, and r are therefore of the 
form 

v = C n e' m43tt + C li e-*'' i4at +e--» mt (C ls cos 1.215f+0 u sin 1.216f), 
<f> = C 2l e' mmt +C 2i <r*- t4it +<r*™ st (C M cos 1.215* + 0 M sin 1.215*), 
r= C sl e' n<stt + 0 82 <r*-""+e-*<>»«(0 88 cos 1.215i + 0 S4 sin 1.2150- 
The particular integrals for any gust may be represented as I v , 
I r , and. their initial values as I m 2^, In, the deriTative of I* being 
I't with the corresponding initial values I 

If, as before (p. 59), we restrict the possible gusts to those of which 
the functional form is different from any of the four functions enter- 
ing into the complementary functions, the particular solutions must, 
on substitution, annihilate the right-hand members of the differential 
equations, and the relations between the constants C</ of integration 
may be determined from the two equations 



Hence, 
and 

Further, 
and 

Finally, 



(O + 0.248)t> +32.17^-1 15.5r = 0, 
0.894« + 0^ + (70.60 + 27.0)r - 0. 

.27240,! + 32.170* - 1 15.50 n = 0, 
.8940 u + 0 0 n + 28.72 0 M = 0, 



0„ = - 8.3260a, 0„ = .2591 0, 



»■ 



- 8.2940 u + 32.170 M - 1 15.5 0 M - 0, 
.8940,, + O0„ - 575.80„ - 0, 

0„=3.7970 2B 0„="!oO58970 JJ . 

■ .O7880 18 + 1.2150 14 +32.170„- 115.50 W «=O, 
■1.2150„,-.O7880 M +32.170 M -115.50 M = O, 
.8940 U +3.920 W +85.740 M -O, 
.8940 M - 85.740 w + 3.920 M = 0, 



and 



0 19 =1O410 M +564.80 M , 
0 U = -564.80„+ 10410,,, 



0 8t =-6.3710 M +lO.560 M , 
0„= - 1O.560„-6.3710„. 
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The solutions therefore, so far as concerns the complementary 
function, are 

4> = <7 J1 e' 02<3 "+ Grf-*"* +c- S2W «7 M cos 1.215*+ O u sin 1.215i), 

v=- 8.326C 31 e- 01 * 3 « + 3.797<7„<r*-*« +<r- SM "[(1041 tf M 

+564.8CJ cos 1.215f+ (- 564.8 #.«+ 1041 C M ) sin 1.2153, 

r =0.2571 O^e- 04 "" +0.005897C a 6-*- Ma{ + e-^f- 6.371 tf M 

+ 10.56OJ cos 1.215f-(10.56J M + 6.371) sin 1.215*]. 

These equations determine the relative magnitudes of the various 
sorts of natural motion. 

The first term is the slowly amplifying divergence, this machine 
being slightly unstable laterally If a side gust is such as to induce 
a lateral velocity of — 8.326 C^, it induces a bank of Q^, an eighth 
as much in radians or seven tunes as much in degrees. It is there- 
fore clear that only very small values of <7 a are admissible for safety. 
The second term, corresponding to the rapidly damped motion, 
shows such rapid damping that it can hardly be of importance, 
except for possible strains on the mechanism, unless 0„ is so large 
that the whole work is inapplicable because of the failure of the 
motions to be small. 

The trigonometric terms show that the oscillation in v will be of 
great amplitude compared with that in <f>, the factor being about 
1200 when d> is in radians or 20 when is in degrees; even the oscilla- 
tion in r will be over 12 times as great as in 6. In other words, the 
machine may have a large oscillatory side-slip or angular velocity 
of jaw without much bank, but for the divergent motion the bank is 
a serious matter for even moderate side-slip. 

The initial conditions <j>=p=v=r=0 give 

0 = <7 M + <7„ + <?„+/*,, 

0 = .02436 ZL - 8.542%, - .3268 G„ + 1.215 + /'*,, 
0 = - 8.326 On +3.797 (7 a + 1041 <L +564.8 U K + /«, 
0 = .2571 <7 n + .005897 0 M - 6.371 C n + 10.56 ff 2i + /„. 

These equations must be solved for the four constants O. 

O a =- .9839/*, -.11487'*, + .0007407 w - .027977™ 
CL=- .000149/*,+ .1170/'*,- .00003427s,- .011637 w 
0~= -.015957*,- .002153/'*, - .000706/ CT + .0396/„, 
Cf t = .01468/*, + .001466/'*, - -0004537/ ro - .0720l7 ro . 

The equations from which the particular solutions are obtained 
are (since Y V =N P = Y r = 0) : 



Ay = (D8 11 -A)v l +L p 6 2 & 1 + (ZAt+iWn 
&4> =7?5i 2 w, + Lpd^ + (X^j, + SM^u 
Ar =/?5 ia y 1 +£jta + (Ms + MJ r i> 



(14) 



422 BEPOBT NATIONAL ADYISOBY COM3UTTEE FOE AEBONAUT1CS. 

or 

Aw = ( - MOD 3 - 95222?' - 346102? + 854)^ + (71342? 

+ 2732)^ - (1 1 25602?* + 1 1047002?)r lf 
= ( - 59.552? - 26. 55)2?^ + ( - 2215025* - 13970Z? 

-34510)^ + (3S952?'+9702?+3062)r 1 
Ar= (-32.842? -280.7)2?^ + (-9030)?! 

+ (- 9922?» - 87242? s - 21032? + 854)r„ 

with 

A = 25922?* + 237802?" + 1800027» + 346102? - 854. 
Motion - in Lateral Gusts. 

We shall take as before the type J(l— e - ") for that of a single 
gust. 

Case 1 . — Side-gust — sharp. v 1 =J{l — e -5 *) . 

2, = J(- 1 + .01473<r«), /„= - .98527/, 

/*=/(- .001028)<r« /*, - - .001028/, 

2'*=/(.00514)<r« /',* = . 00514/, 

I f =/(.002706)<r« / ro =.002706,/, 

0 s1 =-.000384/, = . 0005364/, <?»~ o .000809/, C*= .0002445/. 

The equations of motion are 

lOOOtf/J" = - .384e- M<3 * + .536<r , - M2t - 1.028<r« 

+ e- as8St (.809 cos 1.215« + .2445 sin 1.2150- 

This is all negligibly small. For the same reason certain terms may 
be neglected in v and r. 

v/J*= .OOSe- 0 * 8 * + .002<r* M * - 1 + .01473^ 

+ <T Mm (.98 cos 1.215f-.2022 sin 1.2150, 
100r//= - .01e- 0M6, + .271e-«-«- 3 « s '(.257 cos 1.215* 

+ 1.009 sin 1.2150. 

The effect of the sharp side-gust is to carry the machine sideways 
with it, but not very powerfully at first — much of the air blows 
through the machine — the dominating term at first being 

v=-.2«7«- SM8t sin 1.215*; 

after a few seconds the dominating term is i>= — J, with the very 
slowly growing divergent term effective only after a considerable 
time. There is a slight yawing oscillation, but the extreme angle 
of yaw is only about 0.01/ radians or //2 degrees — the angle being 
computed as 

100*// = J 1 lQQr/J-dt =.4(1- e*™*) + .054 (1 - «r") - .83 16 

+e- !e6 %8316 cos 1.215J + .0122 sin 1.2150- 

The actual sidewise velocity is compounded of v and the amount 
— 115.5* due to the yaw. Hence 



y= j\v- 115.5*)(Zf. 
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For this calculation v and ^ may be simplified to 

v=^/+/e- K ««(cos 1.215f-.2sin 1.215f), 
100*- -.378/-.4/fl- om "+/«-- SM "U32 cos 1.215f); 

and 

y~-.5QJt+l8.5J(e- OIisat -l)-.im 

+/g- SM " (.146 cos 1.215f + .066 sin 1.215*). 

From this it mil be seen that the oscillatory motion is, so far as 
concerns the lateral displacement, of very small amplitude. The 
first two terms which are progressive, are the ones which count. 
Moreover, the displacement is of the same sign as / although the 
side-slip v is of the opposite sign. This apparent contradiction is 
due to the change in orientation the machine moves away from 
the gust owing to the lateral excess wind pressure, but turns into the 
gust owing to the moment of the pressures, and by virtue of the 
great forward velocity-, this turning more than makes up, in the dis- 
placement, for the side-slipping. 

Case 2. — Side gust — mild. v t =J (1 —«—*')• 

/,=/(- 1 + 1.0205<J~"), 7 fQ = .0205J, 

7,=/(.0004043*-"), I* =.0004043 /, 

I\ =J(- .0000809e-«), I'<* = - .0000809/, 

7 r =J(-.001514e~ 2 '), /„= -.001514/, 

O n - - .000331/, tf„ = .00000738/, C„ = - .0000807/, 

C it =. 0001055/. 

It is again seen that there is practically no rolling motion pro- 
duced by the side gust. For v ana r, 

v//= .0027«- otu « - 1 + 1.0205e-» 

+<r-*««(-.0244 cos 1.215f +.1554 sin 1.215i), 
100r// = - .OOSSe-o" 8 " - .1514e-» 

+<r~ SJ «"(.1628 cos 1.215*+ .0672 sin 1.216*). 

(The check «=0, r=0, when *=0, shows that the accuracy has 
been reduced so that the third place is not sure.) The effects of 
the gust are qualitatively as before. The osdflatorv motion is not 
pronounced; the ultimate side-slip velocity is — /; the ultimate dis- 
placement has the same sign as / because the divergent term in 
v — 115.5* is positive. 

Case S. — Side gust — oscillatory. When one examines the records 
made or making at such an observatory as Blue Hill for gustiness 
in the air, no phenomenon is perhaps more striking than the reason- 
ably periodic side switehing of a reasonably steady wind. A south 
wind, for example, may whip back and forth between southsoutheast 
and southsouthwest for hours at a stretch, as Prof. Alexander 
HcAdie has been kind enough to show me on some of his records. 
In the absence of rotary motion, concerning which I am unable to 
find satisfactory data, the simplest way to figure this change in 
direction is as a periodic side gust. A machine goingsouth in such a 
wind would experience an alternating side gust. (The oscillations 
in the head-on velocity of the wind would be relatively very small 
except for actual changes in head-on velocity superimposed upon the 
changes in direction.) It is therefore especially interesting to dis- 
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cuss a periodic side gust — this being the only periodic gust of which 
we can reasonably be said to know anything at all definite. 

Let v=Je i V. We may assume, from our work above that the 
rolling motion will be small and that the side-slip velocity v will not 
be of as much importance in determining the path as the angle f 
coupled with the large forward velocity. The complex value of r is 

(280.7 - 32.84j>i)pWP' 

*" =, 2592p« - 18000^ - 854 + i(34610;p - 23780^) 

If at any one place the period of the complete oscillation is 2t/ti 
with the wind velocity V, the distance traveled by the wind during 
the time of an oscillation in direction is 2vV/n l and this is the dis- 
tance between the nodes of the motion. The time required for this 
machine ( Z7= — 1 15.5) to pass over the distance 2t Vfn is 2t 7/1 15.5n. 
The periodicity of the gust as it appears to the operator of the machine 
will therefore correspond to the value p = 115.5n/V. For instance, 
if F=20 and the tune of an oscillation at one spot were 10 sees, 
so that 71=0.63, the value of p would be about p=3.6, and the 
oscillations would appear to the pilot as taking place about every 
If seconds. A slower oscillation, i. e., a longer periodic time, would 
diminish n and p t — an oscillation at one spot every half minute 
corresponds to a value p= 1.2 on the basis of the assumptions made 
above. 

In considering the values of p which make the amplitude of r large, 
the only hope is to make the term 34610p— 23780j? tolerably small. 
This means o 3 must be about 1.5. For this value, the modulus of r 
is about .03 J and the modulus of the yawing oscillation corresponding 
will be about .025 J. If a wind of 20 foolrseconds is whipping through 
an angle of 45°, the side gust will be only of about 7 toot-seconas 
semi-amplitude and the angle of yaw will be in the neighborhood of 
0.175 radians or 10°. There is nothing to indicate that this would be 
fatal, though it would surely be a nuisance. 

Owing to the fact that the coefficients of i in both numerator and 
denominator are relatively small, the angular velocity I r would be 
about in phase with the gust v lt and hence the angle 7^ would be 
about quartered in phase. If there were periodically an angle of 
10° or 12° between the direction of flight and the relative wind, we 
Bhould find that we were getting into a region where considerable 
rolling and pitching might be induced— for as Hunsaker has shown 
(loc. cat., p. 62) the lateral and longitudinal motions are not strictly 
independent: but as the machine makes the major part of the rela- 
tive wind, the direotions of flight and of the relative wind never 
differ greatly — only some 3° at most in the case under consideration. 

It seems hardly necessary at this time to go into the calculation 
of the actual motion: enough has perhaps been accomplished in 
showing that the oscillation of the direction of the wind induces at 
most a moderate yawing of the machine. The semiamplitude of 
115.5^ would be, if J =7 foot-seconds, about 20 feet; the oenter of 
gravity of the machine would sway back and forth aoross the line of 
night with a total amplitude of 40 feet, until the divergent term 
became effective. 
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Case 4. — Eolling gust. p,= J(l— e - ^). If there were no inter- 
action between v, p, r, the effect on rolling of a rolling gust would be 
figured from the equation 

36. 7Dp+3Up= -314/(1 -«""), 

pfj= -8. 055«r«-««J' e a - 0Mt (l-e- rt )&, 

This means that for any ordinary sharp gust p rapidly acquires the 
value —J, and <f> the value — Jt (radians). It must therefore be 
expected that unless / is very small indeed, the motion will be much 
disturbed^ There will be developed a component of the weight, 
inducing side slipping, and yawing will rapidly develop — the machine 
apparently goes off on a spiral dive. 
We may make the calculations in detail when r = 1. Here 

IJJ= -3. 14-.114<r«, I,o= -3. 25/, 

7 <s /«7'=40.4-l.l<s-*, J*o-39. 3/, 

r*/J=l.l<r<, I'^=1.U, 

Irl 'J '= 10. 58 - . 234«" { , I„ = 10. 35 J". 

tf ai --39. 1/, <?„=. 0027/, tf M =-.219/, tf M =-.163/. 

The equations of motion become 

4>/J= -39. le- t2iset +. 003<r«- 54St +40. 4-1. l<r* 

+(r -^«w(_. 219 cos 1. 215*-. 163 sin 1.215*), 
«//=324e-* 38< +. OlOSe-****-^ 14-. 114er* 

+<!-"«*(_320 cos 1. 2154+49. 2 sin 1. 21K), 
r/J = - 10e Maet + 10. 58 - . 234<r* 

+<r -«6S{(_ 33 L 2l5i+3. 35 sin 1. 215*). 

In the equation for $ the effective terms are 

4>U= -SSie-* 14 **- 1) = -f(nearly), 

and there is a steady divergence in £ to the approximate amount 
— Jt as foreseen. The side-ways velocity v develops more slowly, 
perhaps, but after one seoond amounts to something like 300 J. It 
is clear that / must be very small or the motion becomes disastrous. 

It would be of especial interest to know what sorts of magnitudes 
for / are likely to arise in flight under normal conditions. In so far 
as experience shows that machines are not liable to roll and side- 
slip, it is pretty good evidence that aerial rotary motion with axis 
parallel to the earth is rare and small. 

Case 5. — Yawning gust. r L =J (i—er 1 ). 

IJ J =25.676-*, /„,= 25.67/, 

1*1 J = - .0792 + .154e~', - .075/, 

r*[J= - .154e-', .154/, 

I r [J l-.1235e-S £--.1.12/, 

<7 M =-.006J, C B ~-.00QJ, <7 M =-.0634/, <7 M -.0702/. 
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In this case the motion is 

4>jJ = - .006e- 024Mi - .006e- ww - .0792 + . 154<r* 

+e -JB««(_. 0634 cos 1.215f+.0701 sin 1.215*), 
v/J = + .Obe-* 1 * 3 *' - .023e-*- 5 * ,t + 25.67^ 

+er Ma »(- 26.35 cos 1.215* + 108.9 sin 1.2150, 
r /J= _l_.1235«- < + «-- 3268< (1.145 cos 1.215J + .222 sin 1.2150. 

For moderate values of /, there is nothing serious indicated. The 
coefficients of the divergent terms are small. There can not be much 
roll The most noteworthy phenomenon is the large amount of 
side slip which is fairly rapidly damped out. 

This leaves the rolling gust as the only dangerous type of lateral 
gust. 

The infinitely sharp side gust would produce an initial acceleration 

Yf/J . 

CONSTRAINED AIRPLANES. 

Suppose now that by some automatic steering device the aero- 
plane were constrained to remain pointing in the same direction, 
l. e., so that r=0 identically. The equations of motion become 

(D -Y v )v + (g- YJ))4> = Y*v x + Y p p t + Yp lt 

-L,v+ (kJD-L„)D$=L e v 1 + Lj > p l +Lrr 1 , (15) 

- N v v - N p D<f> = N,v t + ^p^) 1 + N r r 1 + F, 

where Fm is the moment necessary to maintain the constraint. The 
last equation may be regarded as determining F. 

The natural motion of the constrained machine is found from the 
determinant 

A' = « w = 36.72? + 323.1D* + 77.88P +27.15=0. 

This is a cubic equation which has no positive root. 

The negative root is — 8.54. The quadratic factor remaining after 
division by D + 8.54 is 

36.7D»+8.746Z>+3.18 = 0, 

of which the roots are 

Z?=-0.119±0.269i. 

The real part is negative and hence the motion is dynamically stable. 

The introduction of the automatic device has removed the insta- 
bility in the lateral motion. As compared with the complex roots 
in the free motion, these roots indicate a much slower period and a 
considerably smaller damping. 

On the other hand suppose that the constraint had been such as to 
keep the machine level, l. e., 0=0 identically. The equations would 
have been 

(D - Y v )v + ( U- Y r )r= Y v v t + Y p p,+ 3> u 

-1^0-1^-1^+1^,+^ + F, (16) 

- N,v+ QcJD - N r )r=N^ + N p p, + Njr t . 
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The natural motion, would have been determined by 
A"=5 a =70.6l> , +44. 50+109.9=0. 

The roots are 

Z?=-.0.315±0.237t. 

The machine is again Btable. 

It follows that at high speed this Curtiss tractor, which is laterally 
unstable when free, becomes quite stable when constrained either to 
remain on its course or to fly on even keeL 

If stabilizers against rolling and turning were provided, the motion 
would reduce to 

<P-YJv-Y A +Ya l +Yr„ (17) 

and would be stable, Z>= !%=• -0.248. 

It would be a relatively easy matter to discuss the effect of gusts 
of various types on the airplane constrained in various ways; 
two equations are much easier to handle than three. Until some 
definite problem is proposed as important, until some particular con- 
straining device is indicated as likely to be adopted, it may be as well 
not to go into the calculations, which are quite straightforward. 

That a constraint against rolling might be worth while, and would 
indeed be very valuable if rolling gusts were a common thine, is 
suggested by the work done on the free machine (sec 42) where 
gustiness was seen not to be very serious except for the rolling gust. 

DISCUSSION OF METHOD. 

I pointed out in my earlier paper that there were several outs 
about my method of treating gusts. First the gusts must be small. 
If they are not tolerably small, flying would be too difficult — so 
that assumption is not wholly unjustifiable. Second, the calcula- 
tions for determining the mdividual equations of motion and for 
determining formulas for the constants of integration are very 
tedious. Third, the numbers are of such various magnitudes that 
the arithmetical operations which must be carried out cut down the 
accuracy of the work a good deal and indeed, unless great care is 
taken, will lead to illusory or incorrect results. This does not ap- 
pear to be due to any very rapid variation of the true results cal- 
culated from varying data but to the mode of computing. 

To onset these inconveniences we have the satisfactory result 
that once the preliniinary calculations are made, many and varied 
types of gusts may easily be treated, and the further valuable result 
that the actual motion for each case is known so that not only the 
initial motion is determined but the whole extent of the motion. 
This last is necessary for any Just appreciation of the effects of 
periodic gusts and resonance, as has been shown. 

For another method of treating gusts reference may be made to 
a recent paper by Brodetsky and Bryan, "The longitudinal initial 
motion and forced oscillations of a disturbed airplane," Aero- 
nautical Journal, London, 20, 1916, pages 139-166, which has already 
been cited in the text. 
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Much may be said for their method of expansion in series — for 
some problems the work is decidedly simpler than with my method. 
It has been my experience, however, that the application of series 
to the motion of any airplane has its own difficulties and com- 
plicated calculations when the motion is to be followed for any 
reasonable length of time and especially if the machine is defined, 
as 1 have always preferred to regard it as denned, by the actual 
coefficients determined by wind tunnel experiments rather than as 
Bryan's skeleton plane consisting of a main front plane plus tail 
plane — even though the results obtained from such a skeleton may 
be extended to more complicated machines by Bryan's invariant 
method. (See his Stability in Aviation, Chap. VI.) 

The question therefore arises whether there may not be some 
way of abridging the calculations leading to the actual motion of the 
machine. Since finishing my work above, I have received the Pro- 
ceedings of the London Mathematical Society, volumelSj 1917, part 6. 
in which there is an article on "Normal coordinates m dynamical 
systems," by T. J. I'A. Bronrwich, in which he develops a method 
of treating the motions of dynamical systems by means of the theory 
of functions of a complex variable. I wish, in closing, to describe 
the application of Bromwich's work to the problem in hand. 

We have to solve for the longitudinal motion equations of the 
type 

(J) - X u )u - X w w - (XgD + g)9=-P 1 e"i 

Z u u + {D- Z a )w - (Z„ + U)D0 = (18) 
- M K u- Mjo + QcflP- 2ffi)8-P^*, 

where u is a real or complex number, the values we have used being 
o, — r, ±pi. We substutite 

where the integrals are loop integrals in the complex plane and 
£, ti, £ are any functions of X. The results are 

' s^f i-Z^+^-Z u )r,-(Z 9 + U)K]e»ax=P 1 e'" } (20) 
J^jT [_ Mj- M a r, + Qc a W- M q \meVd\=>P t e»'. 
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We next set 

(X-XJf- x m v- (X^+gK =P 1 /(x-/0, 

(k-Z w )n - (Z q + U)H-PJfr-p), (21) 
- MJ- M.v+QJV- 3f a X)r=P I /(X-/0, 
and solve for £, ij, f, finding 

t-P Ai +PAi +P At 

* A(X-/i) ' ^ 

>• — PAi +PA» + P% S ta 
f AQi-n) 

A = 34(X* + 8.49X* + 24.5X 2 + 3.385X + .917) . 

Bromwich shows that, if -with, these values of f , v, i" we take the 
loop integrals (19) around a very large circle, the results for u, w, 
6 -will be the solutions for the motion disturbed from rest at the 
position of equilibrium by the impressed forces P. As he points 
out, this integration is equivalent to the sum of the integrals around 
infinitesimal circles about X=/i and about each of the roofs X of 
A=0, that is, the integral is equal to the sum of the residues of 
ijfi", There is no need to calculate any constants of in- 

tegration. Moreover any of the quantities u, w, 0 can be obtained 
without the others. The numerators in $, % f axe already calculated 
in (20 a, 6, c) of page 59. 

We have, for example, for a head gust u u 

.128X g +1.16X»+3.385X+.917 
S== (X-/*)(X+4.18±2.43i)(X+.0654±.187i) % ' K } 

where the double sign stands for two factors, and « 1 =«7'(1— 
to take a particular case. The residues at each point are merely 
the values of the fraction when one of the factors, the one which 
vanishes at that point, is thrown out of the denominator. In the 
first case for 1 =«** we have as residue of £e\'=£: 
atX=/i=0, 

.917 

~ ( + 4. 18 ± 2.43i) (.0654 ± .187i) = 1 ' 

at X= -4.18 -2.43*, 

■128X !l +1.16X a +3.385X+.917 

( - 4.18 + 2.43i) ( - 4.86i) (- 4.12 - 2.43i ± .187i) ' 

at X=-4.18+2.43i, the conjugate imaginary expression. And 
so on. To treat e~* we should have: 
atX=ji=— 1, 

-.128+1.16-3.385+.917 . 
(3.18 ±2.43i) (.9346 ± .187i) ' 

and so on. 
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As the calculation with, imaginaries involving squares, cubes, 
products, and quotients is by no means simple, it is clear that to get 
the solution for u will be reasonably hard work— much harder than 
to find the particular solutions which for the simple gust involved 
only real numbers. It may be admitted that to work any one gust 
the labor will probably be much less than by my method of determin- 
ing formulas for the constants of integration in terms of the initial 
values of the particular integrals. But as far as I can see, Brom- 
wich's method is of no particular advantage if we desire to calculate 
the effects of a large number of gusts J (1 — e -1 *) of various degrees 
of sharpness both head-on, up, and rotary. When we came to cal- 
culate a periodic gust we found that we were involved in powers and 
products and quotients of complex numbers, and it is probable that 
the work we did in finding the particular integrals was comparable 
with that required for the present analysis. 



In continuation of my previous work in gusts as affecting the 
Curtiss tractor JN2, I have discussed: 

1. Periodic longitvMnal gusts. — It was found that, even in the 
case of best resonance with the slow natural oscillation, the motion 
was not much different from that produced by a simple head-on 
gust until after a considerable time (over 14 seconds) had elapsed. 
The amplitude of the forced oscillation (in up and down motion) 
which ultimately became effective was about Eve times the ampli- 
tude of the gust. This was not regarded as serious because true 
periodicity can rarely be maintained in a head gust and because no 
pilot would wait to let its effect reach such a magnitude. Periodic 
up gusts and rotary gusts were considered as not likely to arise. 

2. General theory of resonance. — It was shown that for aeroplane 

froblems resonance meant different things for different problems, 
t was inferred that resonance was unlikely to be particularly serious 
because in all probability its effect would either be small or would 
take so long to become established that the pilot would check it. 

3. Infinitely sharp ousts. — It was seen that the Bhock to a machine 
was mJLvJ and mZ„J for a head gust, and mX w J and mZ^J for an 
up gust. The serious case is tnZ^J, the vertical Bhock in an up gust 
which was about 4 Jig times the weight, more than twice that found 



serious in a machine where Z w was greater than in the JN2. The 
moral: Keep Z n small, clashes with Hunsaker's conclusion 1 that 
lateral stability is incompatible with high wing loading — but Buch an 
antithesis is common. 3 Reference was made to impulsive gusts. 

4. The effect of {he propeller. — The assumption that a constant 
power instead of a constant thrust was delivered did not very mate- 
rially alter conditions of flight. 

5. Lateral gusts. — The general equations were set up and integrated. 
(a) Single side gusts were shown to produce modern side slipping, 

insignificant roll, and moderate yaw. It was seen that the yaw was 
into the relative wind so that the displacement of the machine in 
space was toward the gust despite the side slipping. 



SUMMARY. 




treated. It would be still more 



l " Dynamical stability of aeroplanes," ■Washington, Smithsonian Mlso. Collect., vol, 62^1818, p. 77. 
> "The production of a laterally stable aeroplane Is attendant with many compromises," Himsakor, p. U 
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Q>) Oscillatory side gusts ■were shown to be a common condition 
of flight, to produce moderate side slipping and yawing, but insig- 
nificant rolling. The path of the center of gravity proved to be 
sinusoidal, so far as the forced oscillation was concerned, and of 
amplitude about two or three times the amplitude of the gust. 

(c) Tawing gusts were found to induce a good deal of side slipping, 
but did not appear to be serious. The roll was very small. 

(d) Rolling gusts were seen to put tie machine into a spiral dive, 
and thus to cause a real danger unless the motion were checked 
promptly by the pilot. 

6. Constrained machines. — A device to keep the aeroplane on its 
course or to prevent rolling made the previously unstable machine 
stable. Such a device might be important to reduce the liability to 
the spiral dive in rolling gusts provided such gusts were common 
phenomena in flying weather. 

7. Other methods of treatment. — The Bryan-Bordetsky method of 
imtinl motions and Bromwich's new method of finding the solution 
for a disturbed state without calculating the constants of integration 
were briefly compared with my system of analysis. 



